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ABSTRACT 

We have measured the correlation functions of the optically selected clusters of galaxies in the Abell and the 
APM catalogs, and of the X-ray clusters in the X-ray-Brightest Abell-type Clusters of galaxies (XBACs) catalog 
and the Brightest Clusters Sample (BCS). The same analysis method and the same method of characterizing the 
resulting correlation functions are applied to all observational samples. We have found that the amplitude of the 
correlation function of the APM clusters is much higher than what has been previously claimed, in particular for 
richer subsamples. The correlation length of the APM clusters with the richness 1Z > 70 (as defined by the APM 
team) is found to be ro = 25.4lg'J /i _1 Mpc. The amplitude of correlation function is about 2.4 times higher than 
that of Croft et al. (1997). The correlation lengths of the Abell clusters with the richness class RC > and 
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1 are measured to be tq — YIA_ ll and 21.0_ 2g h~ Mpc, respectively, which is consistent with our results for 
the APM sample at the similar level of richness. The richness dependence of cluster correlations is found to be 
ro = 0.40d c + 3.2 where d c is the mean intercluster separation. This is identical in slope with the Bahcall & West 
(1992)'s estimate, but is inconsistent with the weak dependence of Croft et al. (1997). The X-ray bright Abell 
clusters in the XBACs catalog and the X-ray selected clusters in the BCS catalog show strong clustering. The 
correlation length of the XBACs clusters with L x > 0.65 xlO 44 /i~ 2 crg s _1 is 30.3±|J /i _1 Mpc, and that of the 
BCS clusters with L x > 0.70 xlO 44 h~ 2 erg s" 1 is 30.2±g;g h^Mpc. The clustering strength of the X-ray clusters 
is much weaker than what is expected from the optical clusters. 
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I. INTRODUCTION 

Rich clusters of galaxies are used to study the large- 
scale structure of the universe on scales larger than 
about 10 r'Mpc (h = H / 100 km s" 1 Mpc" 1 ). To 
measure the clustering strength of clusters the spatial 
two-point correlation function (CF) of the Abell clus- 
ters (Abell 1958; Abell, Corwin, & Olowin 1989) has 
been estimated by many authors (Bahcall & Soneira 
1983; Huchra et al. 1990; Postman, Huchra, & Geller 
1992; Peacock & West 1992; Miller et al. 1999; see 
Bahcall 1988 for a review), and has been found to be 
consistent with the power law 

£cc(r) = (£f (1) 

with the correlation length ro ~ 20 ~ 26 h~ 1 Mpc and 
with the index 7 ~ 2. Recently new catalogs of clus- 
ters have been obtained by automated selection from 
the Edinburgh-Durham Southern Galaxy Catalog and 
the Automatic Plate Measuring (APM) Galaxy Survey 
(Lumsden et al. 1992; Dalton et al. 1997). Redshift 
data of these clusters (Collins et al. 1995; Dalton et al. 
1994b) have been used to estimate the CF (Nichol et al. 
1992; Dalton et al. 1992; Dalton et al. 1994a). The cor- 
relation lengths ro measured for these new cluster sam- 
ples have been reported to be 14 <J ro <J 16 /i _1 Mpc, 
which is much smaller than that of the Abell clusters. 



There are also more recent new catalogs such as North- 
ern Sky Optical Cluster Survey based on Digitized Sec- 
ond Palomar Sky Survey (DPOSS) (Gal et al. 2000) 
and Sloan Digital Sky Survey Cut-and-Enhance (SDSS- 
CE) galaxy cluster catalog (Goto et al. 2002). They 
need, however, either more follow-up confirmations or 
spectroscopic observations before being able to be used 
to study the spatial clustering. 

It has been argued that the counting radius r c = 
1.5 h~ 1 Mpc for the Abell clusters is so large that the 
catalog contains serious projection effects, which cause 
artificial line-of-sight correlations (Sutherland 1988; Ef- 
stathiou et al. 1992). It has also been pointed out that 
the intrinsically subjective nature of the Abell cata- 
log can cause problems in homogeneity and statistical 
completeness (Nichol et al. 1992). However, Bahcall & 
West (1992) have claimed that the discrepancy between 
the Abell and the APM cluster CFs can be explained 
by the richness dependence of cluster correlation ampli- 
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tudes ro = 0Ad c , where d c = n c is the mean inter- 
cluster separation and n c is the mean space density of 
clusters. But Croft et al. (1997) have analyzed richness 
subsamples of the APM clusters and argued that there 
is only a weak dependence of correlation amplitudes on 
the cluster richness. 

It is, therefore, necessary to understand why there 
have been disagreements in the cluster correlation length 
ro and on the relation between ro and d c . In the previ- 
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ous studies individual observational samples have been 
analyzed by different authors who used different meth- 
ods of analysis. In the current situation where the clus- 
ter clustering results from different studies do not agree, 
it is important to reanalyze and compare all preciously 
studied samples in an impartial way. In this paper 
we apply methods of calculating and characterizing the 
CFs to all cluster samples in a same way to remove any 
relative biases. It is hoped that in this objective and 
consistent way of analysis we could find the common 
characteristics of, and intrinsic differences in the spa- 
tial distributions of various cluster samples. 

On the other hand, rich clusters of galaxies often 
have diffuse intracluster gas trapped in their potential 
wells. The thermal X-ray flux from the intracluster 
gas, which is heated to temperature of a few 10 7 K, 
is proportional to the square of the ion density, and 
thus is more confined to the center of the clusters than 
the projected galaxy distribution. Therefore, the X-ray 
selected clusters are expected to have negligible projec- 
tion effects (Briel & Henry 1993; Ebeling et al. 1996). 



Ebeling et al. (1996) have cross-correlated the clusters 
in the ACO catalog (Abell et al. 1989) with the ROSAT 
X-ray sources (Triimper 1993) and created the X-ray- 
Brightest Abell-type Clusters of galaxies (XBACs) cat- 
alog. More recently, Ebeling et al. (1998) have pub- 
lished the Brightest Clusters Sample (BCS) which is a 
true X-ray selected sample from the ROSAT All-Sky 
Survey. It should be interesting to compare the clus- 
tering strengths of these X-ray cluster samples with 
those of optically selected samples to further extend 
our knowledge on the clustering properties of clusters 
of galaxies. 

II. DATA 

(a) The Abell Clusters 

Postman et al. (1992) have published a complete 
redshift catalog of 351 Abell clusters with richness class 
RC > 0, declination S > — 27?5, and tenth-ranked 
galaxy magnitude toiq < 16.5. We restrict our full 



CORRELATION FUNCTIONS OF CLUSTERS OF GALAXIES 



3 



Table 1 

Richness subsamples of the Abell and the APM clusters 



Richness 


N c 


(h 3 Mp C - 3 ) 


ro 

(ft-iMpc) 


7 


Abell 


RC > 
RC > 1 


232 
110 


1.74 x 10~ 5 
0.68 x 10~ 5 


17.411;? 
21.0±*;f 


1 89+ - 15 
1 - oa -0.16 

1 gi+0.34 

-0.36 


APM 


n > so 


218 


2.62 x 10~ 5 


15 4 +0 - 8 

±0 - -0.8 


2 49+ - 17 


■R > 60 

n > 70 

11 > 80 


139 
73 
42 


1.66 x 10" 5 
0.93 x 10~ 5 
0.57 x 10~ 5 


22.3l^ 

25-4l 3 3 ;J 

26.81™ 


1 96+ 016 

2 01+ 29 

2 -o3±8:SI 



sample to 233 clusters with |6| > 40° to avoid serious 
incompleteness at low galactic latitude (Efstathiou et 
al. 1992). In this region with solid angle O = 3.48 sr, 
clusters are nearly uniformly distributed over Galactic 
latitude. Two richness subsamples with RC > and 
RC > 1 are adopted. 

Assuming the Einstcin-dc Sitter universe, we trans- 
form redshifts into comoving distances by using the re- 
lation 

d(z) = 6000(1- l/VTTz) /i^Mpc . (2) 

Since we use the same d(z) relation to all cluster cat- 
alogs, the resulting CFs are not biased with respect 
to one another due to the choice of cosmology. Fur- 
thermore, most other authors have also adopted the 
standard Einstein-de Sitter model in their calculation. 

To estimate the selection function of the survey in 
the radial direction, we first compute the absolute mag- 
nitude M w corresponding to m w of each cluster at a 
rcdshift z from the relation 

m-M = 51og[d(z)(l + z)} + 25 + K(z) . (3) 

The K-correction K(z) is 1.122,2 for the red magnitudes 
of the northern Abell clusters and 4.14z — 0.44z 2 for the 
visual magnitudes of the southern ACO clusters (Ebcl- 
ing et al. 1996). Then the maximum distance D max ^ 
to which each cluster could be observed for the survey 
magnitude limit mio.iim — 16.5 can be calculated. And 
the selection function is given by 

= § E (4) 

i ; max.i 

Using the distribution of clusters smoothed in the red- 
shift space as the radial selection function does not 
make much difference in our results. We limit the 



sample to the distance range from di n = 50 h 1 Mpc 
to d ou t = 500 /i _1 Mpc, and the resulting sample in- 
cludes 232 and 110 clusters for RC > and RC > 1, 
respectively (see Table 1). The region closer than 
din = 50 h~ 1 Mpc is removed to improve the complete- 
ness of the Abell cluster sample and to avoid the possi- 
ble local effects. Large variation of the locations of the 
inner and outer boundaries makes negligible effect on 
our results. 

We estimate the comoving space density of clusters 
in each sample using the method described in Section 
2.3 of Efstathiou et al. (1992). The mean space density 
of the Abell cluster sample with RC > is 1.74 x 
10~ 5 /i 3 Mpc~ 3 while Postman et al. (1992) 's estimate 
is 1.2 x 10~ 5 h, 3 Mpc~ 3 and Efstathiou et al. (1992)'s 
is 1.9 x 10~ 5 /i 3 Mpc~ 3 . The density of our RC > 1 
subsample, n c — 6.79 x 10~ 6 /i 3 Mpc~ 3 , agrees well with 
Bahcall & Soneira (1983) 's 6 x 10~ 6 /i 3 Mpc~ 3 and with 
Efstathiou et al. (1992)'s 7 x 10~ 6 /i 3 Mpc~ 3 . 

Figure 1 shows the distributions of clusters on the 
sky with measured redshifts in each catalog used in our 
CF analysis. 

(b) The APM Clusters 

Dalton et al. (1997) have published a catalog of 
APM clusters with richness range 1Z > 40. This is the 
catalog that Dalton et al. (1994a) have used to mea- 
sure the APM cluster CF, which is consistent with the 
results of Dalton et al. (1992) and Croft et al. (1997). 

The survey region is shown by Maddox et al. (1990)* 
and contains 957 clusters with m x < 19.4 (m x is the ap- 



*The exact APM survey field definition can be found in the 
Stromlo-APM rcdshift catalog (Loveday et al. 1996) served by 
the Astronomical Data Center (http://adc.gsfc.nasa.gov/). 
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parent magnitude of a;-th bright member galaxy where 
x = 1Z/2.1; this is analogous to Abell cluster's m\o). 
The number of clusters in the APM sample shows a 
rather strong dependence on declination (Park & Lee 
1998). In addition to the APM survey limits we have 
restricted our APM samples to the declination range 
of —25° > S > —65° to reduce the potential effects of 
this declination dependence. But this change of decli- 
nation boundaries turned out to make an insignificant 
difference in the CF. 

Since the APM galaxy catalog is complete over the 
magnitude range 17.0 < bj < 20.5 (Maddox et al. 1990) 
and the cluster richness counting slice is [m x — 0.5, 
m x + 1.0], the APM cluster catalog is complete over 
magnitude range 17.5 < ra x < 19.4. This complete 
sample contains 927 clusters. But, in this sample, the 
number of clusters with measured redshifts is only 350. 
For clusters with K > 50 and 17.5 < m x < 19.2, 
however, the completeness becomes 83% and increases 
for higher richness limits. We therefore make four 
magnitude-limited (17.5 < m x < 19.2) richness sub- 
samples with 1Z > 50, 60 (completeness 93%), 70 
(96%), and 80 (100%)(see Table 1). 

To calculate the selection function of the APM clus- 
ters, the bright magnitude limit m Xy i ow — 17 .h should 
be taken into account as well as the faint limit m x , upp — 
19.2. The absolute magnitude M x corresponding to 
m x of each cluster is computed by equation (3), where 
K(z) = 3z is adopted (Dalton et al. 1997). Then the 
selection function of the APM survey is given by 



S(d) 



3 y 



i 



D 3 

max,% 



D 3 . - 

min.i 



(5) 



where D max ^ is the maximum distance to which an 
i-th cluster can be observed by the faint end of the sur- 
vey magnitude m x ^ upp = 19.2, D m i n ^ is the minimum 
distance to which the cluster could be included in the 
sample of m x .i ow — 17.5, and O = 1.02 sr is the solid 
angle of our survey region of the APM sample. Again, 
using the distribution of clusters smoothed in rcdshift 
space as the selection function makes little change in 
our results. The sample is limited to the distance in- 
terval from di n = 50 /i _1 Mpc to d out = 500 /i^ 1 Mpc. 

The mean space densities of our APM 1Z > 70 and 
80 subsamples agree well with those estimated by Croft 
et al. (1997). But the density of our 1Z > 50 sample 
is somewhat lower than 3.4 x 10~ 5 h 3 Mpc~ 3 of Croft 
et al. (1997) because of lower completeness of clus- 
ters with published redshifts. We adopt Croft et al. 
(1997)'s value in the study of the richness dependence 
of correlation amplitudes in § IV. 

(c) The XBACs 

The XBACs catalog (Ebeling et al. 1996) is a com- 
plete, all-sky, X-ray flux-limited (/ x > 5.0 xl0~ 12 erg 
cm~ 2 s- 1 in the 0.1 - 2.4 keV band) sample of 276 
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Fig. 2. — X-ray luminosity L x versus rcdshift z distribution 
of 276 XBACs clusters and 201 BCS clusters. The dashed 
lines represent the L x limit corresponding to the survey flux 
limit; 5.0 X 10 -12 erg cm~ 2 s _1 for the XBACs and 4.4 X 
10~ 12 erg cm~ 2 s _1 for the BCS, respectively. The solid lines 
show the L x limits and distance limits of volume-limited subsam- 
ples. 



Abell clusters of galaxies compiled from the ROSAT 
All-Sky Survey. We constrain our sample to |6| > 30° to 
reduce the effects of the galactic obscuration. After the 
LMC and the SMC regions are removed, the remaining 
survey area covers 6.27 sr of the sky. The complete- 
ness of the XBACs catalog is about 80% due to vari- 
ous sources of incompleteness which are very difficult 
to remove (Ebeling et al. 1996). We have further re- 
moved the X-ray sources contaminated by point sources 
or without measured redshifts from the XBACs catalog. 
The six double clusters in the XBACs catalog are con- 
sidered as single systems which have undergone merg- 
ing. The resulting final sample includes 186 clusters 
with the flux-limit of / x > 5.0 x 10~ 12 erg cm~ 2 s _1 
and with |6| > 30° in the redshift range < z < 0.2. 
This is 91% of clusters in the catalog with the same 
flux, galactic latitude, and redshift limits. 

Assuming the Einstein-de Sitter universe we calcu- 
late the X-ray luminosity L x of clusters from the ob- 
served flux f x and redshift z. Four volume- limited sub- 
samples with L x > 0.24, 0.48, 0.65 and 1.00 xlO 44 hr 2 
erg s _1 and with the corresponding rcdshift limits of 
0.066, 0.092, 0.107, and 0.132 are made (see Fig. 2). 
These samples contain 49, 67, 59, and 51 clusters, 
respectively. Our volume-limited XBACs subsamples 
contain many rich Abell clusters with RC > 1. 

(d) The BCS 

Recently, the ROSAT BCS catalog has been pub- 
lished (Ebeling et al. 1998). Unlike the XBACs catalog 
which is based on the Abell cluster catalog, the BCS 
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Table 2 

Volume-limited subsamples of the XBACs and the BCS clusters 



Lx 

(10 44 ft-^rgs- 1 ) 


N c 


n c 

3 Mpc- 3 ) 


ro 

(ft-iMpc) 


7 


XBACs 


L x > 0.24 
L x > 0.48 
L x > 0.65 
L x > 1.00 


49 
67 
59 
51 


3.53 x 10~ fi 
1.84 x 10~ 6 
1.07 x 10~ 6 
0.53 x 10~ fi 


25.7t 3 3 r 8 
2 5- 2 -4.3 

45 8+ 26 ' 4 


cc+0.51 
z - oo -0.43 

2.48l°;« 

9 91 +0.72 
z - z± -1.27 

1 84+ - 75 


BCS 


L x > 0.24 
L x > 0.48 


33 
49 


3.02 x 10~ fi 
1.70 x 10~ 6 


33.0±g;g 


, 09+O.49 


L x > 0.70 


40 


0.84 x 10~ fi 


30.2±g;g 


9 99+O.94 


L x > 0.90 


34 


0.51 x 10- 6 




2.051J;? 4 



clusters are selected purely by their X-ray properties. 
The BCS contains 201 X-ray brightest clusters of galax- 
ies with measured rcdshifts in the northern hemisphere 
(5 > 0°) and at \b\ > 20°, and is 90 % complete at 
the X-ray flux limit of 4.4 x 10~ 12 erg cm~ 2 s _1 in the 
0.1 ~ 2.4 keV band. Six sources, identified as three 
double clusters in the Abell catalog, are considered as 
six independent X-ray clusters since the sample should 
be based only on the X-ray information. Adding these 
clusters and excluding three clusters with fluxes lower 
than the flux limit, we have 201 X-ray clusters. 

We have used the BCS to make four volume-limited 
subsamples with L x > 0.24, 0.48, 0.70, and 0.90 x 10 44 
h~ 2 erg s- 1 and with the redshift limits of 0.070, 0.099, 
0.118, and 0.134 (see Fig. 2). These subsamples contain 
33, 49, 40, and 34 clusters, respectively. Our high L x - 
limited subsamples of the XBACs and the BCS clusters 
contain very rare objects (see Table 2). 

III. CORRELATION FUNCTIONS 

The correlation function is an estimate of excess 
clustering over the random Poisson distribution. We 
estimate the CFs using the conventional estimator 
DD/DR. We have also used the Hamilton (1993) 's es- 
timator 

, . _ DP ■ RR AN c N r 
Scc{r)- DR2 (A r c _ 1)(7Vr _ x) A . W 

which is less affected by uncertainties in the selection 
function for £ cc < 1 compared to the DD/DR- method. 
We find that the differences between two estimators are 
insignificant in our analysis. Here DD is the number of 
pairs at separation r in the sample with N c clusters, 
RR is the number of pairs in a random sample with 



N r random points, occupying the same survey volume 
and with the same selection function as the real clus- 
ter sample, and DR is the cluster-random pair count. 
We compute the uncertainties of £ cc using the formula 
S£cc = (1 + fcc)/\/DD. Table 3 lists the CF data we 
have calculated. 

(a) Optically Selected Clusters 

Least-square fits of a power law equation (1) to the 
best estimates of CF are presented in Table 1 for the 
richness subsamples of the Abell and the APM clus- 
ters. We determine the 68.3% confidence limits of each 
parameter individually assuming that the distribution 
for the uncertainties of £ cc is Gaussian and different £ cc 
bins are independent of each other, which is a popular, 
but strictly speaking, imprecise assumption. 

Figure 3a shows the CFs of the Abell clusters with 
RC > (circle) and RC > 1 (square). The dashed 
line is Postman et al. (1992)'s power law fit with 
r = 20.6 /i~ 1 Mpc and 7 = 1.86 estimated for the 
Abell RC > sample. The CFs in Figure 3b are those 
of the APM clusters with TZ > 50 (filled circle), K > 60 
(triangle), TZ > 70 (square), and TZ > 80 (open circle). 
The dotted line is a power law with r = 14.3 /i _1 Mpc 
and 7 = 2.05 (TZ > 50) and the dashed line is with 
r = 16.6 fc^Mpc and 7 = 2.10 (TZ > 70) both de- 
termined by Croft et al. (1997). It can be noted that 
the amplitude of CF steadily increases as the richness 
increases. 

Dalton et al. (1992) have compared the richness 
of the ACO clusters with that of the APM clusters 
(hereafter APM92) and found that the mean difference 
in richness is (i?ACO ~ 7^apm92) = 10.1. And Dal- 
ton et al. (1994a) have compared a new version of the 
APM cluster catalog (which is the catalog we use) with 
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Fig. 3. Correlation functions of clusters of galaxies, (a) Magnitude-limited samples of the Abcll clusters with RC > and 1. The 

dashed line is £(r) = (r/20.6 fc^Mpc) -1 - 86 of Postman et al. (1992) for RC > Abell clusters. (6) Magnitude-limited samples of the 
APM clusters with TZ > 50, 60, 70, and 80. The dotted line is £(r) = (r/14.3 h^Mpc)' 205 of TZ > 50 sample of Croft et al. (1997) 
and the dashed line is £(r) = (r/16.6 fc _1 Mpc) -210 of their TZ > 70 sample, (c) Volume-limited samples of the XBACs clusters with 
L x > 0.24, 0.48, and 1.00 xlO 44 h~ 2 erg s" 1 . The dashed line is £(r) = (r/24.6 ft-^pc)- 1 - 80 of Abadi et al. (1998) for clusters with 
L x > 0.54 XlO 44 h- 2 crg s" 1 . (d) Volume-limited samples of the BCS clusters with L x > 0.24, 0.48, and 0.90 XlO 44 /i~ 2 crg s -1 . 



their old version catalog, and found a relation 72-apm = 
1.272-APM92 + 24.8 with a scatter of 16 in 72-apm- We 
then get a rough relation 7£apm = 1.2i?ACO + 12-7. 
Therefore, the counts of the Abell clusters with R = 30 
(RC = 0) and 50 (RC = 1) correspond to the rich- 
nesses of the APM clusters with roughly TZ = 49 and 
73, respectively. In Table 1 we find that the correlation 
lengths of the APM clusters are generally similar to 
those of the Abell clusters at similar levels of richness. 
Contrary to Croft et al. (1997), the APM subsamplcs 
with higher richness limits show clustering as strong as 
or rather stronger than the Abell cluster samples. It 
should be noted that there are large uncertainties both 
in the richness relation and in the richness itself (the 
uncertainties of Raco and 7\Lapm is about 18 and 5, 
respectively; Dalton et al. 1997), and that the Abell 
RC > sample is an incomplete sample with fewer 
poorer clusters. 

Park & Lee (1998) have calculated the CFs of the 
APM clusters using various CF calculation methods 
and sample corrections. They have used different selec- 
tion functions, varied sample boundaries both in angle 
and depth, and applied galactic latitude or declination 
dependence corrections. Even with these variations 



they have consistently obtained CFs with amplitude 
higher than that reported by Dalton et al. (1994a) or 
by Croft et al. (1997), and could not find the sources 
for discrepancy. Dalton (1998) has pointed out that 
our redshift sample contains the clusters with redshifts 
drawn from the Las Campanas Redshift survey or from 
the Edinburgh-Durham-Milano cluster redshift survey 
which cover only parts of the APM Galaxy Survey re- 
gion, and that these clusters could produce an inho- 
mogeneity in the distribution of clusters on the sky. 
Our statistical samples, however, contain very few of 
those clusters, and we have found no effect on the CF 
by them. For example, our magnitude-limited sample 
with TZ > 70 does not contain any of the clusters in 
question. 

(b) X-ray Clusters 

Figure 3c shows the CFs of the XBACs samples with 
L x > 0.24 (triangle), 0.48 (circle), and 1.00 xlO 44 hr 2 
erg s _1 (square). The dashed line represents Abadi et 
al. (1998)'s fit £(r) = (r/24.6 /i^Mpc)" 180 for the 
XBACs clusters with L x > 0.54 xlO 44 h^crg s" 1 . The 
correlation lengths reported by Abadi et al. (1998) 
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Table 3 

Two-point correlation functions of clusters of galaxies 





Abell 










APM 










RC > 




RC > 1 




TZ > 50 




n > 60 




K > 70 




TZ > 80 


r° 


(,cc(r) 


r 


Scc(r) 


r 




r 




r 


«cc(r) 


r 


icc(r) 


3.6 


19 ± 7 


5.3 


12 ± 8 


3.6 


42 ± 10 


3.9 


28 ± 12 


5.3 


21 ± 12 


6.3 


15 ± 14 


5.3 


10 ± 3 


9.5 


4.7 ± 2.2 


5.3 


12 ± 4 


5.8 


12 ± 5 


9.5 


7.5 ± 2.9 


12.5 


7.1 ± 3.1 


7.7 


3.8 ± 1.5 


17.0 


1.8 ± 0.6 


7.7 


4.8 ± 1.2 


8.8 


6.5 ± 1.7 


17.0 


2.4 ± 0.8 


24.7 


0.83 ± 0.59 


11.1 


2.3 ±0.7 


30.5 


0.42 ± 0.20 


11.1 


3.0 ± 0.6 


13.2 


3.2 ± 0.8 


30.5 


0.66 ± 0.24 


48.8 


0.64 ± 0.50 


16.1 


1.5 ± 0.3 






16.1 


0.77 ± 0.26 


19.9 


1.2 ± 0.3 


54.8 


0.19 ± 0.28 






23.4 


0.61 ± 0.16 






23.4 


0.44 ± 0.13 


30.0 


0.61 ± 0.15 










33.9 


0.23 ± 0.08 






33.9 


0.12 ± 0.07 


45.2 


0.13 ± 0.16 











XBACs BCS 



L x > 0.24 b L x > 0.48 L x > 1.00 L x > 0.24 L x > 0.48 L x > 0.! 



r 


€co(r) 


r 


€cc(r) 


r 


£cc(r) 


r 




r 


€cc(r) 


r 


£cc(r) 


6.3 


46 ± 19 


5.3 


57 ± 32 


8.1 


37 ± 27 


5.3 


79 ± 41 


6.3 


52 ± 29 


6.3 
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4.5 ± 2.4 
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15 ± 7 


18.4 
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9.5 


6.8 ± 6.0 


12.5 


5.1 ± 3.4 


12.5 


8.4 ± 9.3 


24.7 


1.5 ± 0.6 


17.0 


1.3 ± 1.2 


41.8 


1.3 ± 0.6 


17.0 


2.6 ± 1.8 


24.7 


0.87 ± 0.73 


24.7 


1.5 ± 1.8 


48.8 


0.19 ± 0.16 


30.5 


1.0 ± 0.5 






30.5 


1.3 ± 0.6 


48.8 


0.32 ± 0.24 


48.8 


0.77 ± 0.59 






54.8 


0.20 ± 0.17 






54.8 


0.44 ± 0.23 











a Separation in units of h 1 Mpc. 

b X-ray luminosity in units of 10 44 /i _2 erg s -1 . 



agree well with our results, except for the highest L x 
subsample with L x > 1.00 xlO 44 h~ 2 erg s _1 . It should 
be considered, however, that the uncertainties in the 
fitting parameters are very large both for their and our 
highest L x subsamples. The CFs of the BCS clusters 
are shown in Figure 3d. 

In Table 2 we list the parameters of CFs of the 
XBACs and the BCS clusters. Our X-ray cluster sub- 
samples have very low number density and consist of 
the rarest clusters which might be very massive. As 
expected, the correlation length is very large and tends 
to increase as the luminosity limit increases. Bor- 
gani, Plionis, & Kolokotronis (1999) calculated the 
CFs from the XBACs and their resulting ro from the 
sample with the same flux-limit as ours, f K > 5.0 x 
10~ 12 erg cm~ 2 s _1 , is 26.0 h~ 1 Mpc, which agrees ex- 
cellently with our result of the lowest L x subsample; 
r = 25.7 h^Mpc. 

IV. RICHNESS DEPENDENCE 

We compare the amplitudes of CFs of the Abell and 
the APM clusters having the same level of richness. 
Richness level of a cluster sample is often represented 

— 1/3 

by the mean separation d c = n c rather than richness 
itself (Szalay & Schramm 1985; Bahcall 1988; Bahcall 
& West 1992; Croft et al. 1997) because the definition 
on richness is different for different catalogs. 



Our results on the richness dependence of the clus- 
ter CF are shown in Figure 4. Filled circles are the 
relations measured from the Abell sample, and filled 
squares are from the APM sample. We find that the 
clustering strength of the APM clusters is consistent 
with that of the Abell clusters, contrary to what Croft 
et al. (1997) have claimed. On the other hand, our re- 
sults agree well with the previous studies on the Abell 
clusters (see Table 1 of Nichol et al. 1994 and other ref- 
erences presented in Fig. 4), but disagree with those of 
the APM team and Nichol et al. (1992). The solid line 
in Figure 4 is a linear fit to the ro versus d c distribution 
found from our optical cluster samples. The relation we 
found is r = 0.40e? c + 3.2. The long and short dashed 
lines are the relations of Bahcall & West (1992) and 
Croft et al. (1997), respectively. Our relation is very 
similar to Bahcall & West (1992) 's law except for a 
constant offset. However, our result strongly disagrees 
with the slope found by Croft et al. (1997). 

The ro versus d c relations of the XBACs and the 
BCS clusters are consistent with each other, but are 
much lower than the relation for optical clusters. The 
slope of the relation is uncertain, but tends to in- 
crease as d c > 80 /i _1 Mpc. The correlation lengths 
calculated by Borgani et al. (1999) from the XBACs 
and by Collins et al. (2000) from the ROSAT-ESO 
Flux-Limited X-ray (REFLEX) galaxy cluster survey 
(Bohringer et al. 2001) are consistent with our re- 
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d c (h" 1 Mpc) 

Fig. 4. Dependence of the correlation length ro on the richness or the mean intercluster separation d c . The solid line is ro = 

0A0d c + 3.2, which is a best fit to our results for the Abcll and the APM samples. The long dashed line is the relation reported by 
Bahcall & West (1992), and the short dashed line is by Croft ct al. (1997). 



suits and their distributions show rather weak depen- 
dence with L x like ours, or look like even randomly 
distributed. As discussed in the following section, we 
think the X-ray luminosity of rich clusters is not a very 
good measure of richness or mass. It should be also 
noted that the uncertainties in r are rather large. 

V. DISCUSSION 

We have measured the cluster-cluster correlation 
function by using four catalogs of optical and X-ray 
clusters. We have applied the same analysis method 
and the same method of characterizing the resulting 
CFs to all observational samples. Our CF for the Abell 
cluster sample is consistent with the previous works by 
Batuski & Burns (1985), Ling, Frenk, & Barrow (1986), 
Retzlaff et al. (1998), and Miller et al. (1999). The cor- 
relation lengths reported by Postman et al. (1992) and 
Bahcall & Soneira (1983) are slightly larger than ours. 
Those estimated by Huchra et al. (1990) and Peacock 
& West (1990) are very close to ours for the Abell clus- 
ters with RC > 1, but are somewhat larger than ours 
in the case of RC > 0. On the other hand, our CF 
for the APM cluster sample has the correlation length 



significantly larger than those reported by Dalton et al. 
(1992, 1994a) and Croft et al. (1997), in particular for 
richer subsamples. We have not been able to find the 
reason for this discrepancy (Park & Lee 1998). 

Sutherland (1988) and Efstathiou et al. (1992) 
have suggested that projection effects of the Abell 
cluster catalog caused an artificial radial correlation, 
which could lead to an overestimation of CF relative 
to the APM sample. They have attempted to cor- 
rect the Abcll sample for the projection effects, and 
reported that the correlation length was reduced to 
r ~ 14 h~ 1 Mpc. However, several other studies have 
pointed out that the uncertainties inherent in this cor- 
rection are very large, and that the correction may have 
been overdone (Nichol et al. 1992). Miller et al. (1999) 
found that the excess radial correlation is not artificial 
but due to two elongated superclusters, Ursa Majoris 
and Corona Borealis. Postman et al. (1992) have esti- 
mated the magnitude of the projection contamination 
and have found that only 17 clusters (5%) in their 351 
Abell clusters may be contaminated. They also have 
measured the effects of these contaminated clusters on 
the CF and have found that the amplitude and slope 
of the CF do not change significantly. Our study indi- 
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cates that the APM team's CFs of the APM clusters 
are underestimated for richer subsamples, and the CF 
of the Abell sample is consistent with that of the APM 
sample at similar richness. 

We have detected a rather strong dependence of the 
clustering amplitude on the richness of clusters. The 
correlation length r versus mean separation d c rela- 
tion we have found from the Abell and the APM sam- 
ples, tq = 0A0d c + 3.2, has a slope equal to that found 
by Bahcall & West (1992), but has a finite correlation 
length at zero separation. We think this is more rea- 
sonable because the correlation length should be finite 
at the continuum limit, namely for the matter field. 

By using the X-ray selected clusters one can greatly 
overcome the projection effects which might have af- 
fected the optical cluster samples selected from the dis- 
tribution of galaxies projected on the sky. Although the 
XBACs clusters are not genuine X-ray selected clusters 
and have incompleteness due to the missing poor ACO 
clusters, they are indeed three-dimensionally bound 
systems and the sampling is nearly free from false de- 
tection or overestimation of richness. The BCS clusters 
are selected purely from their X-ray properties, and 
also have negligible projection effects. However, since 
the number of X-ray clusters is not yet large enough 
to make a sufficiently large volume-limited sample, the 
CFs measured from the XBACs and the BCS clusters 
have large uncertainties compared to those of optical 
clusters. The clustering strength of these X-ray clus- 
ters is in general lower than what is expected from the 
ro versus d c relation of optical clusters. This is proba- 
bly due to the fact that the X-ray luminosity of X-ray 
clusters is not an excellent measure of the cluster rich- 
ness or mass. The X-ray luminosity of a cluster depends 
on the environment and activity as well as its dynami- 
cal mass. In fact, Figure 4 of David, Forman & Jones 
(1999) shows that the clusters very bright in X-ray be- 
long to various richness classes. Gilbank, Bower, & Ca- 
standcr (2001) also found that X-ray selection misses 
some significant rich clusters while optical selection can 
detect all X-ray clusters. However, there is a trend of 
the clustering strength of the X-ray clusters increasing 
as the mean separation increases. Despite the large un- 
certainties, the correlation length of the X-ray clusters 
is higher in amplitude compared to Croft et al. (1997)'s 
expectation. 

We conclude that the Abell and the APM clusters 
do have CFs statistically consistent with each other and 
show a similar richness dependence relation of the clus- 
tering strength. It should be again emphasized that 
our results have been obtained by applying the same 
method of analyzing the observational data and calcu- 
lating the CF. 
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